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We have studied the quasielastic '^He(e, e'p)d reaction in perpendicular coplanar kinematics, with 
the energy and momentum transferred by the electron fixed at 840 MeV and 1502 MeV/c, respec- 
tively. The ='He(e, e'p)d cross section was measured for missing momenta up to 1000 MeV/c, while 
the Atl asymmetry was extracted for missing momenta up to 660 MeV/c. For missing momenta 
up to 150 MeV/c, the measured cross section is described well by calculations that use a variational 
ground-state wave function of the '^He nucleus derived from a potential that includes three-body 
forces. For missing momenta from 150 to 750 MeV/c, strong final-state interaction effects are ob- 
served. Near 1000 MeV/c, the experimental cross section is more than an order of magnitude larger 
than predicted by available theories. The Atl asymmetry displays characteristic features of broken 
factorization, and is described reasonably well by available models. 

PACS numbers: 21.45.-l-v, 25.30.Dh, 27.10+h 



Microscopic calculations make it possible now to cal- 
culate the bound-state and scattering-state wave func- 
tions from Hamiltonian models for processes involving 
three- nucleon systems 111. Thus, using modern (non- 
relativistic) Faddeev and variational |^] techniques 



to solve the three-body problem, one hopes to test 
the ability to predict the structure of three-body sys- 
tems with state-of-the-art realistic NN potentials. The 
quasielastic 3He(e, e'p)d reaction has been used to study 
the single-proton wave function in ^He. In the Plane- 
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Wave Impulse Approximation (PWIA), this reaction 
samples the single-particle momentum distribution in the 
■^He nucleus. However, reaction-dynamics processes such 
as final-state interactions (FSI), two-body currents (me- 
son exchange and isobar), as well as relativity have to 
be taken into account in the data interpretation. Un- 
fortunately, the above mentioned modern computational 
techniques are not yet suffiently developed to reliably de- 
scribe the reaction dynamics at high energies. As such 
they would tremendously benefit fi'om data for guidance 
in their development process. 

High-energy electron beams with high currents and 
100% duty factor at the Thomas Jefferson National Ac- 
celerator Facility (JLab) enable experiments to reach new 
kinematic domains and levels of precision in utilizing the 
(e, e'p) reaction to study nuclear structure and reaction 
dynamics. In this Letter, we address the interplay be- 
tween nuclear structure and reaction dynamics by provid- 
ing an extensive and precise data set that includes cross 
sections and the Atl asymmetry for the ^}ie{e,e'p)d re- 
action in constant quasielastic electron kinematics. This 
data set significantly extends the available data in both 
the transferred four-momentum and the recoil momen- 
tum of the undetected deuteron (missing momentum), 

Prn • 

Measurements were performed using an incident beam 
of 4806 MeV and the two high-resolution spectrometer 
system (HRS) in Hall A of JLab. A detailed description 
of the Hall A instrumentation is available in Elec- 
trons (protons) were detected with the left (right) HRS 
respectively, HRS-L and HRS-R. Scattered electrons were 
detected at a central scattering angle of 16.4° and a cen- 
tral momentum of 3966 MeV/c, corresponding to the 
quasielastic knockout of protons from the ^He nucleus 
with transferred three-momentum \q\ = 1502 MeV/c, 
energy uj — 840 MeV, four-momentum = 1.55 GeV^, 
and Bjorken scaling variable xb — Q"^ / {2u)Alp) — 0.98. 
The range in accepted and oj was ± 0.12 GeV^ and ± 
20 MeV respectively. The ejected proton was detected in 
coincidence with the scattered electron in coplanar kine- 
matics over a range of angles and momenta (see Table 1), 
to measure the pm dependence of the ^He(e,e'p)d cross 
section on both sides of the momentum-transfer direc- 
tion. 

A cooled, 10.3 cm-diameter '^He gas target was used at 
temperature T = 6.3 K and pressures P = 8.30 - 10.9 atm, 
corresponding to densities p = 0.0603 - 0.0724 g/cm^. 
Relative changes in the target density were monitored by 
observing changes in the rate of singles events in the fixed 
HRS-L per unit beam charge passing through the target. 
The target density was determined by measuring the elas- 
tic '^He(e, e) cross section at a beam energy of 644 MeV 
{0e = 30.7°, Q2 = 0.11 GeV^, and normahzing it to the 
cross section derived from a fit to the world data of '^He 
elastic form factors The overall normalization un- 
certainty of the ■^He density is estimated to be 2.9%, ob- 



TABLE I: Central kinematic values of the '^He(e, e'p)d mea- 
surements. Listed are central settings of the hadron (HRS-R) 
spectrometer (momentum Pp, angle Op, and missing momen- 
tum Pm). Negative (positive) pm corresponds to the detected 
proton forward (backward) of q. The electron kinematics were 
fixed, at incident and scattered electron energies of E=4806 
MeV and E'=3966 MeV, respectively, and scattering angle 
of 6*5=16.4° (Q2=i.55 GeV^ |(f|=1502 MeV/c, lo^MQ MeV, 
2;s=0.98). 
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59.83 


425 


1444 


64.76 


550 


1406 


69.80 
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78.28 


1000 


1171 


89.95 



tained as the quadratic sum of the systematic uncertainty 
of our "^He elastic cross section measurement (2.4%), the 
statistical uncertainty (0.5%), the uncertainty in the '^He 
form factors (1.5%), and a 0.5% uncertainty due to pos- 
sible fiuctuations in the target density during the change 
of the beam energy from 4806 to 644 MeV. 

Event triggers were formed by coincident signals from 
scintillator arrays. Particle tracks were reconstructed 
using the HRS vertical drift chambers. The small 7r~ 
background in the HRS-L was rejected using a CO2 Gas 
Cerenkov detector. In the HRS-R, coincident 7r+, ^H, 
and '^H were separated from the protons using the time 
difference between particles detected in the two spec- 
trometers. Most of the accidental coincident events were 
rejected by cuts on the difference between interaction 
points in the target along the beam as reconstructed 
by the two spectrometers, \zh — Zf,\ < 2 cm, where the 
interaction-point resolution was about 8 mm (FWHM), 
and on the E„i (missing energy) spectrum. The £",„ reso- 
lution was 2.4 MeV (FWHM). The remaining accidental 
background was subtracted using the coincidence tim- 
ing between the spectrometers. Events originating in the 
target Al walls were rejected by requiring reconstructed 
events to originate within 3.5 cm from the target center. 
With these cuts, the signal/noise ratio in the most ex- 
treme kinematics, for pm ~ 1 GeV/c, was 50/1 and in 
the worse case, for pm w -600 - -700 MeV/c, about 0.8/1 
- see Fig. 1. 

In the cross-section analysis, a flat acceptance region 
of both HRSs was defined using an R-function cut im- 
posed on the target variables. An R-function is a func- 
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tion whose sign is completely determined by the signs 
of its arguments 0, • Using constructive-geometrical 
properties of R-functions, one can define a complicated 
multidimensional acceptance region as an analytical ex- 
pression, and vary the region's boundaries until the phase 
space is maximized within the flat acceptance region of 
the spectrometers 0. The use of R-functions allowed 
us to double the accepted phase space compared to the 
commonly used rectangular cuts on target variables. 
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FIG. 1: Measured ^He(e,e'p)rf cross section as a function of 
the missing momentum, pm- Also displayed are PWIA and 
full calculations in the diagrammatic approach for two differ- 
ent ground-state wave function. 

The ^He(e, e'p)d cross section was extracted using the 
simulation program MCEEP taking into account the 
effects of internal and external radiation, particle energy 
loss, deviations from monochromaticity of the beam, and 
spectrometer resolutions. For each pm bin, the simulated 
yields were varied by modifying the spectral function 
used in MCEEP to achieve calculated cross sections that 
agreed with the measured ones in both the ^He(e, e'p)d 
Em bin and the adjacent '^He(e, e'p)pn Em bin 9J. Cross 
sections were extracted from the re-weig hted 3He(e, e'p)d 
yield, corrected for radiation, and for contributions from 
'^He(e, e'p)pn to each ^He(e, e'p)d kinematic bin. On av- 
erage, these contributions were about 3%. Within each 
bin, the simulated ^He(e, e'p) cross section was assumed 



to depend on the CTccI prescription of de Forest fo^' 
the off-shell electron-proton cross section. This technique 
allows one to separate the Pm dependence of the reac- 
tion from the rapid dependence on the electron kinemat- 
ics 9]. In addition to the over-all normalization uncer- 
tainty (2.9%, see above), the over-all systematic uncer- 
tainty was 3.4% dominated by uncertainties in the solid 
angle (2.0%), the selection [Em cut) of the two-body 
break-up reaction channel (1.5%) and the knowledge of 
the effective target length via a cut on the interaction 
vertex location (1.4%). 

The extracted ■^He(e, e'p)d cross section is plotted in 
Fig. 1 as a function of Pm- We note that the range of 
Pm measured (resulting in measured cross-section val- 
ues varying over six orders of magnitude), is significantly 
larger than in any other previous measurement. More- 
over, contrary to previous experiments [l3. Iisl H^. our 
measurements over this entire range were performed at 
fixed electron kinematics. 

Also displayed in Fig. 1 are four theoretical curves by 
Laget. The PWIA and full Hannover calculations use the 
Hannover bound-nucleon wave function |lfil | correspond- 
ing to the solution to the three-body Faddeev equation 
with the Paris NN potential and no three-body forces. 
The AV18-HUIX curves are the same PWIA and full cal- 
culations respectively, but with a bound-state nuclear 
wave function derived by a variational technique using 
the Argonne V18 NN potential and the Urbana IX three- 
body force jlS] . All calculations use a diagrammatic ap- 
proach. The kinematics as well as the nucleon and me- 
son propagators are relativistic, and no restricted angular 
(Glauber type) approximation has been made in the var- 
ious loop integrals. Details of the model can be found 
in [l3|. The PWIA curves include only one-body inter- 
actions, while the full calculations include FSI, meson (tt 
and p) exchange and intermediate A formation currents 
as well as three-body (three nucleon tt double scatter- 
ing) amplitudes. The FSI in these calculations follow 
a global parameterization of the NN scattering ampli- 
tude, obtained from experiments in LANL, SATURNE 
and COSY (18|. On the scale of Fig. 1, the differences 
between the calculations using the two ground-state wave 
functions are very small. By far, FSI constitute the major 
difference between the full and PWIA calculations. Me- 
son exchange and intermediate A current contributions 
are generally small (up to 20-25%), and the three-body 
contributions are negligible 0. 

Three regions of Pm can be discerned in Fig. 1. For 
\pm\ below ~ 150 MeV/c, roughly within the Fermi mo- 
mentum, the deuteron can be viewed as only marginally 
involved in the interaction. Hence, the data are expected 
to be dominated by the single-proton characteristics of 
the '^He wave function. As can be observed, both the 
PWIA and full curves describe the data quite well, and 
the difference between them is rather small - see also 
Fig. 2 for details. For \pm\ between 150 and 750 MeV/c, 
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well above the Fermi momentum, the cross section is ex- 
pected to be dominated by the dynamics of the reaction. 
Indeed, very large contributions from dynamical effects 
are observed. While the full calculations describe the 
data very well, the PWIA curves over-predict the data 
by up to a factor of 2 for 150 < < 300 MeV/c and 
under-predict them by up to an order of magnitude for 
300 < < 750 MeV/c. The differences between the 
two PWIA and two full-calculation curves are very much 
dominated by FSI. At xb=^, the on-shell rescattering of 
the fast nucleon on a nucleon at rest is preferred and the 
contribution of FSI is maximal. Because the NN scat- 
tering amplitude is almost purely absorptive in the JLab 
energy range, the corresponding FSI amplitude interferes 
destructively with the PWIA amplitude below, and con- 
structively above pm « 300 MeV/c We note the 
difference in cross sections in this region for negative and 
positive pm, and it is discussed below. For p^ larger than 
750 MeV/c, the calculations gradually deviate from the 
experimental data: at 1000 MeV/c, they grossly under- 
predict the measured cross section by more than an order 
of magnitude. Whether it is a consequence of the trun- 
cation of the diagrammatic expansion or a signature of 
other degrees of freedom is an open question. 
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FIG. 2: Same data as in Fig. 1 for low pm only, but shown 
as a ratio to the full calculations using the ground-state wave 
function (gswf) generated from the AV18 NN potential and 
the Urbana IX three-nucleon force. Also shown are the ra- 
tios to this calculation of the full calculations that use the 
Hannover gswf, as well as of the two corresponding PWIA 
curves. 

The sensitivity of the data to the details of the wave 
function at low \pm\ is shown in Fig. 2. In order to en- 
hance the details. Fig. 2 displays the low \pm\ subset of 
the data from Fig. 1 as a ratio to the full calculations 
using the AV18 NN potential and the Urbana IX three- 
nucleon force. Also displayed are the ratios to the same 
calculation of the full Hannover ground-state wave func- 
tion and the two corresponding PWIA curves. As already 
noted, in the low \pm\ region, we expect reaction effects 
such as FSI and two-body currents to be relatively small 



as compared to higher | (Fig.^ , and hence the data to 
be more sensitive to the details of the calculated ground- 
state wave functions than to the uncertainty in describ- 
ing reaction dynamics. As can be seen in the figure, the 
curves produced by this model are mainly sensitive to 
the details of the bound-nucleon wave function. We note 
that for Pm below 50 MeV/c, the calculations are purely 
co-planar perpendicular kinematics whereas experimen- 
tally, because of the large it is difficult to avoid con- 
taminations with parallel and out-of- plane components. 
For \pm\ > 50 MeV/c, we observe that the curve that 
best agree with the data is the full AV18-I-UIX. We sug- 
gest that this better agreement with the data is related 
to the fact that the wave function generated from the 
AV18-fUIX potentials reproduces the correct ^He bind- 
ing energy (by construction), while the Hannover wave 
function that does not include three-body forces under- 
binds the ^He by ^ 0.7 MeV. 

The Atl asymmetry was extracted for < \pm\ < 
660 MeV/c according to 



TL 



(1) 



where (T+ and a- are coplanar '^He(e, e'p)d cross sections 
measured at positive and negative missing momentum 
respectively. The Atl observable downplays the signifi- 
cance of the ground-state wave function, by virtue of the 
ratio involved in its definition |l9l | and there exist indica- 
tions that it is sensitive to relativistic effects and to 
mechanisms that break the simple factorization scheme 
of PWIA cross sections . 

Figure O displays the extracted A-tl data with the 
PWIA and full calculations using the two ground-state 
wave functions described above. The difference in the 
two ground-state wave functions has a very small effect 
in the full calculations. In contrast to the PWIA calcula- 
tions, the measured Atl displays a structure characteris- 
tic of broken factorization [21| : the oscillating pattern of 
Atl comes directly from the interference between differ- 
ent reaction amplitudes. Both full calculations describe 
the data reasonably well by displaying similar structure. 
Such structure in Atl was previously observed in the 
quasielastic removal of p-shell protons in the ^^0{e,e'p) 
reaction [i^l, and was well reproduced by relativistic 
Distorted- Wave Impulse Approximation calculations by 
Udias et al. In that case broken factorization was 

attributed to dynamical relativistic effects, enhancement 
of the lower components of the Dirac spinors. However, 
these effects are marginal in our experiment because of 
the low nuclear density of '^He [2^. Rather, in our case 
the factorization is broken by the strong interference be- 
tween the PWIA and re-scattering amplitudes . 

In summary, we measured the ■^He(e, e'p)d cross sec- 
tions and Atl asymmetry at = 1.55 GeV^ and 
xb = 0.98. For \pm\ below 150 MeV/c the data are 
mostly sensitive to the details of the bound-state wave 
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FIG. 3: The measured Atl asymmetry. The curves are the 
same four calculations used in Figs. 1 and 2; by definition, 
the two PWIA curves are indistinguishable. 



function. The best agreement is observed with calcula- 
tions using a '^He ground-state wave function generated 
from the Argonne V18 NN potential and the Urbana IX 
three-nucleon force, which also better reproduces the '^He 
binding energy. For \p„i\ from 150 to 750 MeV/c, strong 
FSI effects are observed as quenching (enhancement) of 
the cross section below (above) \pm \ of about 300 MeV/c. 
For missing momenta from 750 to 1000 MeV/c, the mea- 
sured ^He(e, e'p)d cross sections are increasingly larger 
(more than an order of magnitude at 1000 MeV/c) than 
predicted by available theories. Whether it is a conse- 
quence of the truncation of the diagrammatic expansion 
or a signature of the existence of exotic effects is an open 
question. The measured Atl displays strong structure 
characteristic of broken factorization due to interference 
between the PWIA and re-scattering amplitudes. Cal- 
culations using a diagrammatic method well describe all 
observables up to l»mj_= 750 MeV/c. Other calculations 
of this reaction |3j US that have recently become 
available similarly interpret the data. 
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